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<^ ■ We present the concept of selective readout for broadband resonant mass gravitational 

D ■ 

I ' wave detectors. This detection scheme is capable of specifically selecting the signal from the 

contributions of the vibrational modes sensitive to the gravitational waves, and efficiently 
rejecting the contribution from non gravitationally sensitive modes. Moreover this readout, 

> ■ 

■ applied to a dual detector, is capable to give an effective reduction of the back-action noise 



o 



within the frequency band of interest. The overall effect is a significant enhancement in the 



predicted sensitivity, evaluated at the standard quantum limit for a dual torus detector. A 
CO , molybdenum detector, 1 m in diameter and equipped with a wide area selective readout, 



would reach spectral strain sensitivities ~ 2 x 10 ^■^ Hz between 2-6 kHz. 
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{^JQ, The concept of sensitive and broadband acoustic mass gravitational wave (GW) detectors has 

been recently presented with reference to a dual sphere detector . Such a novel class of detectors 
may be of great interest, as they would be sensitive in a broad frequency interval of few kHz 
in the kHz range, were GW signals from fully relativistic stellar sources are expected jj]. A 
practical implementation of such a detector could suffer from: i) unwanted non-GW active resonant 
modes in the frequency region of interest, ii) additional thermal noise, due to the low frequency 
contribution of the high frequency non GW sensitive modes P|. These difficulties may be overcome 
by introducing a geometrically selective readout, capable of specifically selecting the contribution 
to the signal from all the GW sensitive modes. We have found also an additional bonus: the scheme 
is capable to give an effective reduction of the back-action noise within the bandwidth. The overall 
effect is a significant enhancement in the sensitivity. 

In a "dual" detector one would measure the differential displacement, driven by the GW, of 
the nearly faced surfaces of two concentric massive bodies, mechanically resonating at different 
frequencies P|. In this scheme the centers of mass of the two bodies coincide and then remain 
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mutually at rest while the masses resonate. The differential displacement is then only due to the 
internal vibrational modes, while the centers of mass provide the rest frame for the measurement. 
This design allows the use of wide bandwidth (non resonant) readouts, evolution in concept and 

n 

in technology of the resonant readouts used for bar detectors jj| . 

A basic dual detector can be represented as a simple one dimensional system [Fig. H^a)] , where a 
force Fe, acting on two different mechanical resonators, is evaluated by a differential measurement 
of their positions xi and X2- In the frequency region between the two resonant frequencies, Fg 
drives the slow resonator above its resonance Us and the fast one below its resonance I'f. The 
responses of the two resonators are then out of phase by tt radians and therefore the differential 
motion results in a signal enhancement over the single oscillator responses [Fig. ^b)], as shown 
by the transfer function Hp^ = Xd/Fg [Fig. ^c)]. 

Such a scheme in addition leads to a reduction of the back-action noise, just within the detector 
bandwidth of interest. In fact the system response to the back-action force, that it is coherently 
applied with opposite direction on the two masses, is nearly exactly in phase and the consequent 
differential displacement is highly depressed at a frequency /**, as shown in the transfer function 
Hba = Xd/Fha (Fig. ^). We measure x^ with a displacement amplifier, described by its additive 
displacement and back-action force noises, with white power spectra Sxxi^) = Sxx and Sff{uj) = 
Sff, so that the total displacement noise is Sxx + \Hba{^)\'^ S f f . The noise power spectrum on the 
measurement of F^, due to the amplifier system, is then: 

SfA^) = {Sx. + \HM\^Sff)/\HF^{u,)\^ . (1) 

If we take as reference an operation at the Standard Quantum Limit (SQL), we may consider 
SxxSff = and the noise figure can be optimized by adjusting the ratio Sxx/ Sff. In a wide 
bandwidth detection strategy Sxx and Sff must be balanced to give the lowest noise within the 
bandwidth. In doing so we profit by the subtraction effect in the back-action noise transfer function 
i/fea, and obtain a dip at the frequency /**, as shown in Fig. ^d). We finally notice that, to fully 
exploit the back-action reduction features, /** should be placed, by a proper choice of the system 
parameters, amid the oscillator frequencies. 

In the case of a three-dimensional body, the dynamics of elastic deformations is given as the 
superposition of the dynamics of an almost infinite number of normal modes of vibration 5]. An 
obvious way to preserve the convenient features for signal and back action noise outlined above, is 
to bring the real system to be as close as possible to the idealized two modes system. In fact when 
only the first quadrupolar mode can be considered for each body, the response to a GW of such 
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FIG. 1: (a) Onc-dimcnsional "dual" detector: the same force Fe is measured by the relative displacement 
Xd of two resonators, (b) Transfer functions of the slow resonator (continuous line) and of the fast resonator 
(dashed line), (c) Dual detector transfer functions: signal Hp^ — Xd/F^ (continuous line), back-action 
Hba — ^dlFba (dashcd hnc). (d) Wideband optimized noise. 

a system can again be described by the simple one-dimensional model. This can be accomplished 
with a novel selective readout we propose here, capable of rejecting a large number of normal modes 
on the basis of their symmetry. This geometrically based mode selection senses the surface position 
of a body on specific regions, so that the related deformations can be combined with a weight /sign 
properly chosen to optimize the total response to normal modes of quadrupolar symmetry. Such a 
strategy to select specific vibrational modes and to reject a class of unwanted modes is conceptually 
different from the strategy now employed in GW acoustic detectors. Resonant bar detectors a| 

normal modes excited by the GWs 
by the use of resonant displacement readouts coupled to the modes. With a proper choice of the 
read-out surfaces, the resonant scheme is not sensitive to the thermal noise of out of resonance 
modes and gives an efficient frequency based mode selection. However this feature necessarily 
limits the detector bandwidth, due to the thermal noise contribution of the resonant transducer, 
as in every resonant readout scheme. By contrast the geometrically based mode selection selects 
gravitationally sensitive normal modes by means of their geometrical characteristics, and shows a 
significant rejection of non quadrupolar modes without affecting the detector bandwidth. For this 
reason it can be effectively applied to a "dual" wideband detector. 

Without loss of generality, we apply these new concepts to a dual torus detector [Fig. [2fa)], a 
convenient geometry where the advantages of the proposed scheme can be fully exploited. In the 



FIG. 2: (a) The two concentric torii, Tg, Ti are made of materials of density pe, pi and have the same 
height. The inner torus may also have null internal radius and be reduced to a cylinder. The relative distance 
between the two bodies is measured in 4 regions (in black), each of area St^ for the whole torus height, 
(b) Section of the detector showing the signal enhancement obtained when a GW signal drives the external 
torus above resonance and the internal torus below resonance. The difference Xd = xi -\- — X2 — X4, 
proportional to the GW strength, is not dependent on a number of non GW sensitive modes. 

dual torus we average the differential displacement in 4 distinct areas j;i..4 [Fig.|2{b)] and combine 
them to obtain X^i = xi + — X2 — X/j^. The detector displays its maximal sensitivity when the 
GW propagates along the z axis, the symmetry axis of the system. The corresponding force does 
not depend on z and the system response can be well described by plane strain solutions, where 
the displacements are functions of x and y only and the displacement along z vanishes. All the 
components of the internal stress are also independent by z. In this case the displacement normal 
modes of a single torus, eigenfunctions of the free body dynamic equations, are functions of the 
kind: 



where the functions fa,n,g a,n are linear combinations of Bessel functions of the coordinate r, with 
coefficients given by the boundary conditions. The integer a represents the angular symmetry of the 
mode, while n identifies the mode order within the angular family. The orthogonal displacement 
fields w^„,w^„ represent the same radial distribution of the deformation, mutually rotated by 
for this reason they share the same eigenvalue Ua^n, called resonant frequency of the mode. 
Any displacement u may be written as linear superposition of these basis functions: 



s, a, n 

where the time dependent coefficients are determined by the force acting on the body and s = +, x . 



^^tni^) = fa,n COS (06*) V + g a,n sin (a 6*) 
Wa,n(r) = -fa,n Sm{a9)ir + g a,n COS {a 9) iQ 




(2) 
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A + polarized GW propagating along the z-axis applies on the mass of density p the force 
density Fgy,{t) Gg^(r): 

Fgw{t) = \ph{t) 

Ggy,{Y) = r cos(2 6l)ir -r sin(26')i0. (3) 

For symmetry reasons this force can only excite w+ quadrupolar (a = 2) modes. A "weight 
function" approach to the problem will give the mathematical framework to study the selective 
read-out, which implements the difference Xd = xi + — X2 — [Fig.[2Ib)]. If Ug and Uj are the 

n 

displacements of the torus Tg and Tj 0\ , we define the observable physical quantity of the system 
as: 

Xd{t)= [ [ue(r,t)+u,(r,t)] ■ ^^dV, (4) 

where the "selective" measurement strategy and detection scheme is implemented by the weight 
function P4(r) = PJ (r) P^{0) i^, where: 

Pa = 5{r - Te) - 5{r - Ti) 

1 4 
r?i=0 n=0 

and 0(x) represents the unit step function. The normalization is P/v = St, area of one sam- 
pling region [Fig. |2{a)]. Here the term gives the requested displacement difference rg — rj 
[Fig-Htb)], while reduces the angular integral over 4 distinct regions, 2a wide, centered around 
= 0, ^,7r, ^ [Fig- Eta)]. A value a = 0.3 rad is assumed to perform the following calculations. 
For comparison we consider a non selective transducer system, which senses the displacement over 
a single area, centered for example at = 0. Its weight function Pi has angular component 

= Q[e + a) - @{e - a) + Q{e + a-2Ti) - @{e - a- 27r), while the radial dependence and 
the normalization remain the same as P4. 

When we evaluate the observable Eq. (jlj), by using the expansion Eq. ^ and the weight 
functions P4 or Pi, each w^^ mode contribution depends on its angular symmetry a. As shown 
in Fig. El^b), in both cases the modes contribution oscillates and rapidly decreases, due to the 
averaging over the area St- But in the case of P4 only the symmetry values a = 2 -|- 4A:, with k 
non negative integer, give non null contributions to our observable Eq. (0J). To summarize, the 
quadrupolar modes family [k = 0) contribution is essentially preserved, many mode families are 
rejected {a ^ 2 + 4k, k > 0) and the residual higher order families {a = 2 + 4:k, k > 0) give a 
reduced contribution. 
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We point out that the normal modes w^, proportional to s'm(aO), and excited by x-polarized 
gravitational waves, are rejected by P4(r), for every value of a. When a second transducer system, 
identical but rotated by 7r/4, is employed to detect these modes, any z-axis propagating GW 
can be fully characterized in terms of intensity and polarization. 

The detector sensitivity can be evaluated by the transfer function Tx^ = Xri.(co)/F(uj), which 
gives the observable Xd induced by a generic driving force density G(r) We call 

(w^ Wa.n) and (v* zua,n) the normal modes and eigenfrequencies of the torus Tg and Tj, 
while the loss angles {(p, ip), inversely proportional to the mode quality factor Q, describe the 
dissipation. We have: 

>p /G(r)-w^,„(r)^/w^,Jr)-^dy / G(r) • v^,„(r)rfy Jvg Jr) ■ ^Mrfy 

where P(r) can be P4(r) or Pi(r). If F{t) and G(r) are given by Eq. 0, by Eq. (jB)) we calculate 
the system transfer function to a GW, defined as Hg,^{uj) = X(i{u)) /h{u)). In the simple case 
Pe = Pi = Pi we have: 

while the more general case is straightforward. The read-out back-action force is applied in the 
sensing areas with the same intensity but opposite sign for the two bodies, so that its spatial 
density is given by The corresponding system transfer function is then: 

UA^) = TxA^)\G{r)=nr)/PN ■ 

In close analogy with Eq. the transfer functions Eqs. [T] and |H1 give the detector sensitivity 
to GW as: Shh{^) = {Sxx + \Hi,a{io)\'^ S f f) /\Hgw{uj)\'^ . In Fig. |2tc) we compare the sensitivities of 
a detector, evaluated for the selective readout P4 and for the single area read-out Pi. 

The main effect of the selective readout is to cancel out both thermal and back-action noise 
contributions due to the non quadrupolar modes, so that a flat response is obtained in the full 
bandwidth of few kHz and the back action reduction feature is exploited. The detector dimensions 
are properly chosen to profit by the GW sensitivity of the second order quadrupolar mode of the 
outer torus. We also notice that the use of large read-out sensing area highly reduces the thermal 







. In fact, as the shorter wavelength 



noise due to the cumulative effect of all the normal modes 
modes are averaged out, a very good transfer function convergence may be obtained by adding less 
than 100 modes. 
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FIG. 3: (a) Angular dependence of the weight function P4 (continuous line) and of the displacement induced 
by a w+ mode (dashed line) with symmetry a — 2: the mode contribution is proportional to the integral of 
the product of these two functions, (b) Normalized angular contribution of normal modes as a function of a 
evaluated for the weight functions P4 (hollow symbols) and Pi (crosses), (c) Predicted sensitivity of a Mo 
dual torus detector (same as Fig. 0)1 with the selective read-out P4 (continuous line) and with the standard 
read-out Pi (dotted line). 
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FIG. 4: Predicted spectral strain SQL sensitivities of two different dual detector configurations. The pre- 
dicted SQL sensitivities of two advanced interferometric detectors are also shown (dotted lines): LCGT and 
a narrow band design of Advanced LIGO. Continuous line: Mo dual detector, inner cylinder radius 0.25 
m, outer torus internal-external radius 0.26-0.47 m, height 2.35 m, weight 4.8 + 11.6 tons, fundamental 
quadrupolar modes 1012 Hz and 5189 Hz, amplifier noise S"^^ = 5 x 10^^^ mVHz, Q/T > 2 x 10^ K'^. 
Dashed line: SiC detector, inner cylinder radius 0.82 m, outer torus internal-external radius 0.83-1.44 m, 
height 3 m, weight 20.5 - 41.7 tons, 5*:^:,; = 3 x lO^^^^ mVHz, Q/T > 2 x 10^ K'K 

Finite element analysis demonstrates the selective read-out rejection capabilities also for many 
classes of 3-dimensional vibrational modes. The sensitivity enhancement over the standard readout 
scheme is then not limited to our plane strain approximation, and the dual torus configuration could 
be evolved in a complete detector. A practical readout configuration is a series of 4 capacitive 
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transducers, gradiometrically connected and sensed by a single SQUID amplifier. This implements 
the selective scheme P4 and the consequent back-action reduction. The needed sensitivity could 
be reached as recent progress show that the SQUID amplifiers are approaching the quantum limit, 
now also in the necessary strong coupling configuration 1 1211 . and that the electric polarization field 
can be increased up to the material intrinsic limitations (l3|. 

To evaluate the limits of our design, we can evaluate the sensitivity at the SQL of some practical 
configurations of detector material and geometry. As usual a low dissipation material is required 
to reduce the effect of the thermal noise. Molybdenum represents an interesting choice, as it shows 
high cross-section for GWs and its mechanical dissipation was investigated at low temperature 
giving Q/T > 2 x 10^ K^^ for acoustic normal modes ^J]. In Figs. El and |1] is shown the SQL sen- 
sitivity of a Mo detector with dimensions within the present technological production capabilities. 
In Fig. lUwe also show the SQL of a detector made of SiC, a ceramic material currently used to 
produce large mirrors or structures with mechanical and thermal properties of interest here 
but not jet characterized in terms of low temperature mechanical dissipation. 

The selective readout scheme applied to the dual concept allows the design of detectors tailored 
for relatively high frequency GW and with very few spurious modes within their wide bandwidth. 
These features could make the dual torus detector complementary to advanced interferometric 
detectors, as shown in Fig. by the comparison with the expected sensitivities of LCGT J^] and 
one of the possible setting of Advanced LIGO in narrow band operation 
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